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AUmCT 

An  •pproKlanc loo  for  cospuclng  cb*  b«*y  ^to4  cbnrncCMrtMU*  la 
eh*  Cl/C/1  qu«u«,  MUMCcd  bjr  D.Ctom  nad  CJI.lnnrta  la  chair  hook 
‘'ruadaaantaU  of  (>M«alag  Thaary”,  la  ahoim  to  ha  lachlac  la  accuracy, 
and  la  charafora  9«lca  ondaairahla  for  practical  t*aa. 

KET  yptos 

Cl/C/1  quaua,  buay  parted,  fhaaa  Type  dlacrlbaclon,  fhaaa  Typa  ianawal 
Rrocaaa,  coapucac tonal  prohabtllty. 
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l.latgcdttctlont 

la  a racMt  book  en  queualat  cbaory  by  Croaa  aa4  HarrUfll,  aa 
approalaaca  aacbod  la  aastaatad  for  obtaiatac  cba  baay  parlod  cbaractarl#' 
tlca  of  cha  Cl/C/1  quaua.  Tbroutb  a atapla  caaputatlaaal  aaaaipla,  wa 
daaoaatrata  that  cbls  aacbod  caa  ba  aavaraly  laaccaraca.  It  ta  poaalbla 
Co  coaacrucc  (alapla)  aaaaplaa  tdiara  cba  approalaac ioa  la  arbUrarlly 
bad.  Our  dlacuaaloa  wlli  ibou  cbac  cba  approalaacloa  could  raaulc  la 
a groaa  uadar-asclaacloa  of  cba  duraclea  of  cba  buay  partod  aad  la 
cbarafora  uadaalrabla  for  uaa  la  pracclcal  atcuacloaa. 

Wa  will  draw  our  aaaapla  froa  cba  claaa  of  CI/N/l  quauaa  wlcb 
lacar-arrlval  claaa  followlat  a gaearaltiad  Crlaag  dlacrlbucloa.  Tbaaa 
quauaa  baloag  Co  cba  claaa  of  FH/H/1  quauaa,  ^Icb  la  Cure  ara  a a«b> 
claaa  of  cha  M/C/1  quauaa,  aaalytad  la  dacail  by  kaaaauaal  |S|.  la 
Sactloo  2,  wa  provide  a brlaf  dlacuaaloa  of  Phaaa  Typa  dlacrlbutloea  aad 
Phaaa  Typa  ranawal  procaaaaa  aad  follow  cbla  la  SacCioa  3 with  a 
discuaaion  of  cba  FN/H/1  quaua  wlcb  parcicular  rafaraoca  co  lea  buay 
pariod.  la  SacClon  4 wa  dlacuaa  cha  approalaaca  aacbod  auggaacad  by 
Croaa  and  Harris  obcalnlng  cha  spacial  fora  of  chair  approxlaaca 
foraula  for  cha  aaan  langch  of  cba  busy  period.  Finally  la  SacCloa  S, 
wa  present  our  coapucational  exaaple  and  aaka  soae  further  raaarfcs 
on  Che  approxlaacion. 

2.  FH-Dlstrlbutions  and  FM-Ranawal  Frocassaa: 

A phase  typa  probability  distribution  (FH-'disCrlbution),  introduced 
by  M.F.Neut8(2l,  is  obtained  as  the  distribution  of  the  tiae  till 


•bsorpeton  in  « ftalc*  Nsrluw  eliatA  «rtch  m •bsorbUkC  stat*  late  iilitjck 
•baorption  la  cartala.  U*  iImII  briafly  ravlaw  aucb  aodala. 

Cooaldar  a coaclauaus-ctaa  Markov  cbala  wlcb  staco  apaco 

Ut ,o,  oai)  for  vbteb  cbo  otarM  1 a ara  traaalaat  aad 

(iH>l)  is  abaorbtng.  Xc  la  aaauaad  (hat  a(ar(U«  la  aay  craaolaac  acata 
(ha  chain  gats  abaorbad  la  (a^l)  with  probability  oaa.  Xha  laflal- 
taalaal  ganarator  of  auch  a Markov  chain  than  haa  tha  fora 
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%#hara  T la  an  n»n  noo-alngular  aatrla  with  and  V|j^*  Alao 

T®>  0 la  an  n-vactor  aatlafylag  Ta  ♦ T®  • 5,  vhara  a*»(l, 1). 

Ua  aaauaa  that  wa  ara  alao  glvao  an  initial  probability  vactor  (a.  0) 
for  tha  Harkov  chain. 

It  la  now  anally  saan  that  for  tha  abova  Markov  chain  tha  tlaa 
till  abaorptloo  in  tha  atata  (wfl)  haa  c.d.f. 


F(a)  - 1-0  aapCTi)  a,  a jO. 


Any  probability  diatributioa  F(.)  ao  obtalnad  will  ba  callad  a Fhaaa 
Typa  Dlatrlbutloo.  and  tha  pair  (o,  T)  la  callad  a rapraaaotatloo.  Tha 
acan  of  F(,)  la  glvan  by  X ^ ■ -oft  ^a. 

In  [3]  Neuta  considered  Che  renewal  process  obtalnad  by  restarting 
Che  Markov  chain  Q Instantaneously  after  each  absorption  (renewal)  by 

performing  a multinomial  trial  with  probabilities  o^  and  outcomes  1(...( 


•nd  c«ll«d  le  • fha—  Typ»  t<wi^wi»l  froci 


— (PM-t— wwl  ProcMc).  I%i« 
rMulc*  lA  AAotbAf  coAtiAueuA-clAA  Itorkov  chAla  wlch  •tat*  cmc* 

I I a ) aad  lat  lalCMlaal 

Q*  • T ♦ T®  A®, 

tfhar*  A®«dlat(a.( .a  ) aad  T®  • (T^ Harkov 

chala  daacribaa  cha  ‘>haaa‘*  of  cha  PH-raaaaal  ^ocaaa.  la  12|  it  liaa 
baao  ahoMi  that  cha  rapraaaotattea  (a,  T)  of  P(.)  oay  ha  ao  choaaa  chat 
cha  aacria  Q*  la  trraduclbla.  M»  ahall  aaauaa  chat  chla  ta  ladaad  tht 

caaa. 

For  cha  Plh-raoawal  procaaa  daacrlbad  ahova,  lot  Hfc)  doaoca  tha 
nuabar  of  raaaaala  la  (0,c|  aad  J(t)  daooca  tha  phaaa  at  CLoa  c^  aad 
daf loa 

Flj(v,  C)  • P(R(t)*  V,  J(t)-J  I J(0>*1|,  v»  0,  I 2 • •* 

la  (3)  It  la  ahovn  chat  cha  gaaoracing  fuaccioo  of  cha  aatricaa 
P(v,  t)  • t))  ia  glvan  bp 

P(s,  t)  - ^«''p(v,  t)  » a*p{CP«-tT®A®)t)  , 0* 


Also  by  dif faranciacing, 

H(t)  . |_  Hz,  0|^,.  - J,  n «*'’  o'*  * 

where  the  (1»  j)-th  entry  of  H(t)  is 


(2.1) 


- E(  H(t)  1 


■ a 


I JC0)-1) 
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Not*  CluiC 

■(e)  • aH(c)«  • a r — Q T®  (t.t) 

k-l  k! 

is  ch«  •Jip«cc*4  ntaib«r  of  rooovolo  ta  (0»  !)■ 

3.  Tho  W/H/i  Quooot 

Uo  nov  coooldor  • CI/H/1  quooo  !■  whleli  orrtvolo  occur  occor41«t  eo 
• Phoso  Typo  rooouol  procoos.  Socli  • aodol  wttl  bo  doaotoitf  by  Ht/H/l,  Mo 
lot  It  donoto  cbo  oorvlco  roco. 

Lott  log  X(c)  ond  J(c)  doooco  roopocclooly  cbo  guouo  longcb  (l.o.*  cbo 
nuabor  of  cuocouoro  lo  cbo  oyocoa)  ood  cbo  pbooo  of  cbo  orrlool  proeooo 
oc  clao  c-o,  ic  io  oooily  oooo  cboc  ((X(c),  J(c))t  c>0)  lo  o Hortov  cbola 
trlch  ococo  spoco  (0,1, )>(l,....,a)  ood  lof iaitooUail  goaorotor 
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The  inf Inltesiaal  generator  Q Is  of  a fora  studied  by  V.Wallace  (7)  under 
the  naae  quasl-blrth-and-death  procesooo.  A queueing  aodel  with 


infinitesiaal  generator  of  the  saao  fora  has  also  been  studied  by  Meuts  (A) 


For  Ch«  Narkov  chain  Q ahova  w*  lat  dnaota 

cha  probabllley  that,  atartlat  U tha  atata  (!♦!,  j),  1^,  tha  flrat 
alalt  to  tha  aat  of  atataa  i • ((l.h)t  oecara  no  latar  t^  « ana 

that  at  tha  apech  ot  aach  a flrat  paaaasa  tha  phaaa  of  tha  apatan  la  j*. 

Ut 

Cj^.Ca)  • 7 a"**  aCj^.U)  , ta  a 

and  lat  C(a)  aanota  tha  n^o  aatain  of  aotrlaa  C^jiCa). 

Thaoraa  1. I 

1)  Tha  natrlx  C(a)  aatlaflaa  tha  aquation 

C(a)  • (al  - (T-vl))*'  ul  ♦ (al  - (T-aUrVA®  C*(a)  (1.2) 

for  all  a>0. 

11)  If  tha  quaua  la  atabla,  l.a..  If  o«l/li<l,  than  0>€(0«>)  la  atochaatlc. 

Ill)  C la  tha  ainiaat  aolutlon  In  cha  aat  of  aub** atochaatlc  antrlcaa 
to  cha  aquation 

C - (ul  - T)*^  lul  ♦ T®A®C*|  (1.1) 

la)  Tha  aatrix  C can  ba  coaputod  uainf  tha  racurrooca  fornulaa 


Cq  -0 


(1.*) 


Proof: 

Equation  (3.2)  follows  fron  a standard  first  passage  arguaent  bjr  considering 
Che  first  cine  Che  queue  length  goes  either  dcwn  or  up.  That  the  queue 


r 


f 


Is  stsbls  Iff  p>A/u<l  i«  t«sU-lu>owi.  ths  proofs  of  cJio  oclMr  eUlas 
In  ehs  thoorsB  srs  sUillsr  to  choso  in  llsuts(4l  sad  banco  onltcod. 

RMisrkt  In  IS)  le  Is  shown  chse  cba  nstrls  C Is  serirtlp  pesltiva  and  thus 
has  a ualqua  Inwarianc  prebabtllcp  vaccor  f sadsCirlot  if>g»  Aa>l. 

Thaorsn  3.S: 


Than 

ji*  - -(l-C«C)t(T  ♦ T®A®)  ♦ (T®A®  - pl)C)~^a.  (J.A) 

• • 

whars  C is  cha  o«n  nstrlx  sach  of  ahosa  rows  Is  g. 

Proof t 

Dlf fsrsnclacint  <3>2)  with  raspacc  to  s and  lattint  s*0,  wa  gat 

l(T-pl)+T®A®-H®A®C)£*  • -a.  (3.7) 

Sines  1-CMi  is  non-slogular  and  slnca  |•la(T-1ll)0*T®A®C^*0•  «>a  hava  that 

I(T-pl)+T®A®4T®A®Cj(l-C4C)  • (TfT®A®)  ♦ (T-iila2T®A®)C 

Using  this  in  (3.7)  wa  gat  (3.6)  by  noting  that 
(T+T®A®)C  - (Tt>T®A®)^  - Og  • 0 

Rotark:  dote  that  v*  is  tha  expactad  duration  of  a busy  pariod  «diich 
starts  with  1 custonar  and  in  phasa  i. 

Thaoren  3.8; 

gp*  “ l/|v(l-p))  (3.9) 

Proof ; 

This  fomula  is  established  for  the  nore  genet Jil  M/G/1  queue  in  (5). 
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R— ark;  Etjuaelon  (3.9)  providea  a powerful  coaputatlooal  check  on  Che 
accuracy  of  nuaerical  coapucat looa . 

Altorltha  to  coapute  the  aean  leotth  of  the  buay  period; 

Step  1;  Coaiputa  C uains  (3.4) 

Step  2:  Coapute  Che  Invariant  prohabillty  vector  f of  C 
Step  3:  Conpute  ji  ualnt  (3.6) 

Step  4;  Verify  (3.9) 

a 

Step  it  Conpute  the  aean  lencth  of  the  buey  period  E(BP)*a|i 


4.  The  Croea-Harrle  Approxlnatloo; 

In  their  book  |1], Croee  and  Harris  obtain  a “nice,  neat"  expression 
for  the  expected  length  of  the  busy  period  In  a CI/C/I  queue.  By  treating 
the  arrival  process  ss  chough  It  wore  nenoryless  and  nlalcking  Che  derivation 
of  the  ncan  length  of  the  busy  period  In  an  N/C/1  queue,  they  "show" 
that  the  expected  length  of  the  busy  period  In  a CI/C/1  queue  is 
approxlnacely  equal  to 

E(BP)  • E(S)/(1  - 7 a(c)  dB(t))  (4.1) 

0 


where  S denotes  the  service  tine  and  B(.)  its  c.d.f.,  and  where  n(t) 
Is  the  expected  nunber  of  arrivals  in  (0,  t]. 

Note  Chat  for  Che  PH/H/1  queue 

,o  .“Pt  , 


7-(t)  dB(t)  . « 7 ? ^ I pe  dt 

0 0 k-1 


by  (2.2) 


• p o(l-p  Q ) T , 


whence  (4.1)  nay  be  written  as 

E(BP)  * l/(p  - 0(1  - p“V)‘^  1®1 


(4.2) 
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The  existence  of  the  Inverse  In  (4.2)  follows  from  the  fact  that 

A A 

Q Is  a stable  matrix,  l.e. , Re  6j^<0  for  every  eigenvalue  6^  of  Q . 

5.  The  Computational  Example; 

For  our  example  we  considered  a queue  with  a generalized  Erlang 
distribution  for  Inter-arrlval  times  which,  as  a PH-dlstrlbutlon,  has 
the  representation 
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The  service  time  distribution  was  assumed  to  be  exponential  with 
parameter  p > 0.483.  It  may  be  verified  that  for  this  queue  p - 0.900. 

In  Implementing  the  algorithm  at  the  end  of  Section  3 to  compute  the 
(exact)  value  of  the  expected  length  of  the  busy  period,  we  terminated 
the  Iterations  In  (3.4)  using  the  criterion 

l<M<n  < I”'®- 

Then  a linear  extrapolation,  based  on  the  last  two  iterates,  was  used  to 
make  the  computed  matrix  G stochastic.  The  Invariant  probability  vector 
of  G was  computed  using  Wachter's  method  [6].  All  these  were  programmed 
In  APL,  and  a copy  of  the  program  may  be  obtained  by  writing  to  either 
of  the  authors.  Me  sunmarlze  the  results  of  our  computations  below: 


-t- 


- (lA. 55679,  20.58772,  29.53159,  30.85796. 

£ - ( 0.27996,  0.54192,  0.05679.  0.05938, 

- 20.73940 
p-^(l-p)-^  - 20.73940 

a p*  - 14.55679 

The  approximation  to  the  expected  length  of  the  busy  period  wee 
computed  using  (4.2).  The  computed  value  of  the  “approx laaC ion**  for 
E(BP)  turned  out  to  be  5.95040  which,  when  compered  Co  Che  true  value 
14.55679,  is  grossly  inaccurate. 

A close  examination  of  the  model  here  reveals  why  the  “approxlMClon" 
fares  so  poorly.  The  assumption  of  a aMiMryless  property  for  the  arrival 
process  in  the  example  amounts  to  resetting  the  phase  of  the  PH-renewal 
process  at  1 immediately  after  every  departure.  This,  for  our  model, 
has  the  effect  of  “slowing  down"  the  arrivals  limaediately  after  each 
departure  and  leads  to  under-estiiaating  the  duration  of  (and  the  ntaabcr 
served  during)  a busy  period.  What  one  needs  to  recognixc  is  that  the 
resetting  of  the  phase  to  1 is  repeated  several  times  during  a busy  period 
and  not  just  at  the  end  of  the  first  service  completion  alone.  As  p gets 
larger  and  closer  to  1,  one  expects  a busy  period  to  consist  of  a larger 
niimber  of  service  completions  making  the  approximation  progressively 
worse  as  p Increases.  Further  by  considering  a general  matrix  T 


32.24840) 

0.06194) 


of  the  form 


10- 


T • 


-« 


-a  a 


-b  b 


-b  b 


•b 


for  daflnlnt  cha  ganaralliad  Erlaag  dlacrlbucloo  aed  bjr  aakljig  *a'  aaallar 
and  'b*  Urgar.  ona  can  conatmcC  axaaplaa  whara  cba  approataatloo  ia  t#oraa 
for  qnauaa  vlch  avao  cha  aaaa  valua  for  o. 
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Gl/G/I  queue,  busy  period,  Phase  Type  distribution 
Phase  Type  Renewal  Process,  computational  probability. 

A,  


20.  fConttnum  on  rorotoo  old*  I#  nocooooPT  onrf  IrfonMfy  hf  MocM 

An  approximation  for  computing  the  busy  period  characteristics  in 
the  GI/G/1  queue-~ASuggested  by  D.  Gross  and  C.M.  Harris  in  their  book 
"Fundamentals  of ^Queueing  Theory", *ls  shown  to  be  lacking  In  accuracy, 
and  is  therefore  quite  undesirable  for  practical  use. 
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